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We show that the injection of polymer chains into nanochannels becomes easier as the channel 
becomes narrower. This counter intuitive result arises because of a decrease in the diffusive time 
scale of the chains with increasing confinement. The results are obtained by extending the de Gennes 
blob model of confined polymers, and confirmed by hybrid molecular dynamics - lattice-Boltzmann 
simulations. 



The dynamics of confined polymer chains in solution is 
of great importance to many biological processes. These 
range from the permeation of macromoleculcs across the 
nuclear pore complexly] , to their packaging and ejection 
from viral capsids 0, 0L and their motion in the crowded 
cellular environment 

A problem that has captured the attention of the 
physics community in recent years is the translocation 
of polymer chains through narrow pores [f|, an ubiqui- 
tous process in cellular biology that has enormous po- 
tential in developing controlled applications such as the 
separation of macromolecules in designed microfluidic de- 
vices [6 1 or the selective recognition and sequencing of 
DNA [2 1 . While much effort has been devoted to charac- 
terise the transient passage through very short pores @- 
11 1, for long channels [l2[ questions such as whether the 
chain translocates at all, and how this process is trig- 
gered, have been explored to a much lower degree. 

From a theoretical perspective, it is now well accepted 
that the injection of chains into long channels involves the 
surpassing of a free energy barrier [^, 0] , which for flow- 
driven chains is set by the competition of the driving hy- 
drodynamic drag and the cntropic pressure resisting the 
deformation of the chain as it is squeezed into the con- 
striction. Within the de Gennes blob model [HI , Sakaue 
et al. [l4[ showed that the barrier is overcome above 
a critical flow rate, which, remarkably, is independent of 
the length of the chain and the geometrical features of 
the nanochanncl. 

In this Letter we highlight the importance of hitherto 
unknown scaling aspects of the barrier-crossing dynamics 
by exploring the Brownian motion of the chain within the 
channel. We extend the de Gennes blob model beyond 
the current static description, and include the effect of 
the attempt frequency of the chain to cross the free en- 
ergy barrier. Surprisingly, we find that the threshold flow 
rate decreases with decreasing diameter of the nanochan- 
nel, making translocation more probable across narrower 
constrictions. This unexpected result is due to a decrease 
in the diffusion timescale of the blob- like chains as the di- 
ameter of the channel is decreased, therefore increasing 



the attempt frequency of the chain to cross the free en- 
ergy barrier. To confirm this prediction, we have carried 
out hybrid molecular dynamics - lattice-Boltzmann sim- 
ulations (Til , flij , finding a very good agreement with our 
theory. 

This novel and counter intuitive feature of the translo- 
cation process has relevance in nano- and biotechnolog- 
ical applications. While we focus on the flow-driven 
translocation of the chains, the enhancement of translo- 
cation in narrower pores is independent of the driving 
force, and is thus expected to hold for other systems in- 
volving confined polymer chains in solution, for example, 
for electrophoretic jY?! and diffusiophoretic [HI driving. 

Blob model for polymer translocation. - Exact analytic 
or fully resolved numerical approaches to the dynamics 
of a confined polymer chain in solution are, in general, 
not feasible. The chain itself contains a large number 
of atoms and the complexity of the problem is increased 
even further by the need to include the hydrodynamic 
interactions mediated by the fluid. 

To make progress, coarse graining is needed, where 
the number of degrees of freedom is drastically reduced 
in a way that preserves the relevant physics in a given 
problem. Our approach here is to extend the de Gennes 
blob model to study the flow-driven translocation of the 
chains. 

We start by considering a linear polymer composed of 
N monomers joined by links of size 6, which lies con- 
strained by a channel of diameter D > b. We focus in 
the limit of small pores, where the chain is deformed 
from its equilibrium radius of gyration, R ~ bN 3 / 5 3> D, 
adopting an elongated shape of length L. The blob model 
assumes that the chain can be regarded as a sequence of 
M blobs, each of size £ (see Fig. [I]). It then follows that 
the extension of the chain inside the channel is L ~ £M. 
Within each blob, the effect of confinement is assumed 
to be unimportant. Therefore, the size of the blob scales 
with the number of monomers in a blob, P, following the 
Flory prediction, £ ~ 6P 3 / 5 . Using the conservation of 
the number of monomers, N = MP, one can obtain a 
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FIG. 1. (Color online) Blob picture for the penetration of a 
linear polymer chain into a nanochannel driven by a conver- 
gent fluid flow. When pushed into a nanochannel of diameter 
D, the polymer chain deforms into an elongated configura- 
tion composed of a series of space-filling blobs of size £ ~ D. 
As the penetration length, x, increases, the chain can reach 
the position of the free energy barrier, x* , where hydrody- 
namic driving takes over the resisting entropic pressure, thus 
triggering the translocation of the polymer. 



scaling law for the length of the chain: 
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This result is consistent with the Flory free energy of the 
confined chain, 
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where T is the temperature and k B is Boltzmann's con- 
stant. Eq. ([2]) describes the competition between the 
entropic elasticity of the chain conformation, given by 
the first term, and the internal collisions that give rise to 
excluded volume effects, represented by the second term. 
The length of the chain is set by the balance between 
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both contributions, giving the scaling L ~ D (£) ' N, 
which matches the blob model result, Eq. ([TJ, for space- 
filling chains, i.e. where £ ~ D. 

Turning our attention to the translocation process, we 
consider a polymer which is brought to the entrance of 
the constriction by an underlying fluid flow and partially 
enters it up to a distance shown in Fig. [TJ If the free 
portion of the chain is not significantly deformed by the 
fluid pressure, the only competing mechanisms governing 
the motion of the polymer are the entropic cost of confine- 
ment, which tends to increase the free energy of the chain, 
AF, and the drag exerted by the driving fluid, which 
tends to reduce it. Assuming an entropic contribution of 
k B T per confined blob, and that the fluid drag acting on 
each blob scales as ijJ^/D 2 , it follows that AF depends 
on x, according to AF(x) ~ k B T (§) - ^ (§) 2 01; 
here J is the flow rate (measured in units of volume per 
unit time) and 77 is the dynamic viscosity of the fluid. It 
is clear from the interplay between the two terms in AF 
that there is a free energy barrier, AF* ~ (k B T) z /(77J), 



which is overcome once the chain has reached a critical 
penetration length x* ~ ^j-D. 



Sakaue et al. [14[ considered the case where the 



strength of the barrier is of the order of thermal fluctua- 
tions, i.e., AF* ~ k B T. This leads to x* ~ D, showing 
that the translocation is triggered when a single blob has 
moved into the channel. It follows from this assumption 
that the barrier is surmounted above a critical flow rate 
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The striking feature of this leading-order scaling predic- 
tion is that J c depends on neither N nor D, as long as 

l>x* HUHg. 

However, this first approximation disregards the fluc- 
tuating motion of the chain in the constriction, which is 
responsible for a crucial dependence of the translocation 
on the diameter of the nanochannel. Even if J < J c , 
the chain can still overcome the barrier due to thermal 
fluctuations. This is easily seen by considering the prob- 
ability of chain translocation, P(J) — KT m < 1, where r m 
is the observation time and n = kjj exp {— AF*/k B T} is 
the translocation rate over the barrier. The transition 
rate contains the exponential factor indicative of the free 
energy barrier, and the attempt frequency, Kd, which re- 
flects the fluctuation of the position of the chain within 
the channel. This additional contribution to P(J) is par- 
ticularly important close to the penetration threshold, 
where the forces acting on the polymer effectively bal- 
ance each other out leaving the blobs in the chain free to 
undergo a Brownian motion over length scales compara- 
ble to their own size. The mean square displacement of 
the blobs evolves in time, t, according to (x) 2 ~ D blob i, 
where the diffusion coefficient obeys Einstein's relation, 
F>biob ~ k B T/(rj^). In this limit the attempt frequency 
follows by taking the reciprocal of the timescale of self- 
diffusion of one blob, kq — -Dbiob/C 2 - Recalling that for 
a space-filling chain £ ~ D, we get 
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To demonstrate the relevance of this result let us define 
a threshold current, J* < J c , via the relation P(J*) = 
P* , where P* < 1 is an arbitrary value of the probability 
of translocation. After substitution of kd we obtain 
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where X 3 = k B TT m /(r]P*). We thus find a logarithmic 
correction to the leading-order scaling prediction which 
gives the dependence of the translocation rate on the 
observation time, e.g., the experimental timescale, as well 
as on the channel size, D. According to Eq. ([5]), J* is 
only weakly sensitive to D when the observation time is 
long enough. In this limit, the simple scaling structure 
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of Eq. ((3]) is recovered. However, if the observation time 
becomes comparable to the blob time scale ~ k^ 1 , then 
J* increases sharply with D. 

Molecular dynamics /lattice- Boltzmann simulations.- 
In order to confirm the predicted dependence of the 
threshold flow rate on the nanochannel thickness, we have 
carried out numerical simulations of confined chains in 
narrow channels, in equilibrium and in flow-driven con- 
figurations, using a hybrid molecular dynamics - lattice- 
Boltzmann algorithm. 

We consider a bead-spring model for the polymer. The 
chain is coarse-grained into N solid beads of radius a 
joined by N — 1 masslcss springs. The dynamics of the 
chain is given in inertial form, 

mr i = -^V ii f/ + F i! (6) 

where is the position vector of the i-th bead, m is 
its mass and the double dot symbol stands for double 
differentiation with respect to time. The first term in 
the right-hand side contains the elastic force between 
adjacent beads in the chain, which are modelled using 
a Hookcan potential with elastic constant k and rest 
length 6, and the excluded volume interactions between 
all beads in the chain which follow from a DLVO po- 
tential, C/dlvo = U exp{-Krm(rij - 2a)}/{r ij - 2a), 
where m is the distance between beads i and j, kdh 
is the Debye-Hiickel screening length and Uq is an am- 
plitude. The same potential is used to model repulsive 
interactions between the beads and solid walls. 

The remaining term, F, = F\ + F^, accounts for the 
random and viscous forces, FJ and F^, that originate 
from thermal fluctuations and viscous stresses within the 
fluid. The random force, FJ, has zero mean and a time- 
correlation matrix satisfying the fluctuation-dissipation 
relation, (FJ(i)FJ(t')) = 2k B T(5(t - t')l, where C is the 
friction coefficient of a bead with the fluid, S(t) is Dirac's 
delta function and 1 is the identity matrix. The viscous 
force, F* 1 , is modelled using Stokes law, which compares 
the velocity of the bead, r,, and the velocity of the fluid 
at the position of the bead, v(rj); Fj- 1 = — £(fj — v(rj)). 
Here £ = Gnnr^ is the Stokes drag coefficient, where r^ 
is the hydrodynamic radius of the beads. 

The dynamics of the solvent is modelled using the 
fluctuating lattice-Boltzmann algorithm for a Newtonian 
fluid flBl fig , which we solve simultaneously to Eq. ([6]) . 
This allows us to compute the friction force, F' 1 , on 
all beads, which is exerted back onto the fluid to en- 
sure global momentum conservation. The bead mobility 
~ l/Qifqa is recovered by including a lattice correction 
to TV 

Confined chains in equilibrium.- We first compare the 
equilibrium properties of polymers under full confine- 
ment to the blob-model prediction. We consider a chain 
in a periodic simulation box. In lattice units, the di- 
mensions of the box are 240 x D x D along the x, y 
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FIG. 2. (Color online) (a) Length of a confined chain L, as 
a function of the number of beads in the chain, N, for two 
channels of different width, D. For large N and small D, sim- 
ulation results (symbols) approach the blob-model prediction, 
characterised by a scaling exponent v = 1. (b) Length of a 
chain as a function of the channel width. Symbols correpond 
to simulation results showing a good agreement with the ex- 
ponent of —2/3 predicted by the blob model. Solid lines in 
both panels are a linear fit to the data. Error bars are smaller 
than the symbol size. 

and z directions, respectively. Also in lattice units, chain 
parameters are fixed to N = 128, a = 0.25, b = 1.0, 
k = 'SOkftT, and m = 1 7(4a 3 /3)p. DLVO parameters are 
fixed to Uq = k^T and kdh = 80, which ensures a repul- 
sive force between beads comparable to fc^T at distances 
between beads ~ 0.5a, preventing chain crossings. The 
viscosity, density and temperature of the fluid bath are 
chosen as rj — 6, p = 36 and k^T — 0.1, while we take 
the hydrodynamic radius as = 0.32. 

The initial configuration is prepared by placing a 128- 
bead chain (R « 8) in a wide channel (D = 16); we subse- 
quently apply a small forcing along the negative x—, y— 
and 2— directions and allow the chain to relax to a 'com- 
pressed' state which is then used as the initial condition 
in our simulations, reducing the likelihood of trapping 
in metastable states. Equilibration is then carried out 
for several values of D and N over 1.25 x 10 4 timesteps, 
during which we record the length of the chain, L, by 
measuring the difference between the maximum and min- 
imum x— coordinates accessed by the polymer at a given 
time. Data is averaged over fifty realisations of the noise. 

Following previous numerical studies [2l| - [25| . we ex- 
pect a scaling behaviour L ~ N v , where the ^-exponent 
should cross over from the 3D value v w 0.59 to the ID 
value v = 1 as the confinement becomes stronger. Our 
numerical simulations confirm the scaling prediction for 
the confined regime, as can be seen in Fig. [DJa), where 
the upper curve, corresponding to D = 2, already shows 
a linear increase for moderate values of N. For D = 8 we 
observe a slightly smaller exponent of v = 0.89 ± 0.01; 
such proximity to the fully developed ID value, even for 
weakly confined chains (D ~ R w 8), indicates that 
the crossover occurs rather quickly, both as a function 
of D and N. This assertion is further supported by look- 
ing at the dependence of the length of the chain on the 
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FIG. 3. (Color online) (a) Probability of translocation as a 
function of the flow rate for different pore sizes, (b) Threshold 
flow rate, J* , as a function of the channel width. The thresh- 
old flow rate increases with increasing width in agreement 
with the blob model (see text). Error bars are comparable to 
the size of the symbol. 



channel width for relatively long chains (N = 128), which 
shows the expected power-law behaviour with an expo- 
nent — 0.68±0.04, very close to the blob model prediction, 
L ~ D- 2 / 3 (see Fig.^b)). 

Forced translocation.- To study the translocation pro- 
cess, we equilibrate a 128-bead chain for 1 x 10 4 timesteps 
on the cis side of a two-compartment duct of dimen- 
sions 240 x 14 x 10 separated by a channel of dimensions 
24 x D x D. We disregard the process of finding the 
pore, and carry out the equilibration stage while keeping 
the first bead of the polymer tethered to the channel en- 
trance according to the same elastic potential as used for 
the beads. We subsequently release the chain and exert 
a constant body force, /, on the fluid in the x-direction, 
thus creating a convergent flow that drives the chain into 
the channel. Simulations are carried out for a further 
1 x 10 5 timesteps, which allow for the record of failed 
(chain ejected back to the cis side) or successful (chain 
released to the trans side) events. 

Fig. E^a) shows the probability of translocation of the 
chain, P, as a function of the imposed flow rate, J, 
for varying D. Data are in units of a reference flow 
rate, Jo ~ fo, where the reference body force is set to 
/o = 1 x 10~ 5 in simulation units. Due to the increase 
in the drag with J, the probability of translocation in- 
creases with increasing flow rate, crossing over from a 
non-translocation regime (P — S- 0) to a full translocation 
regime (P — > 1). 

While the observed increase of P as a function of J 
holds for all channel widths considered, the curves in 
Fig. Eta) shift to the right as D is increased. This means 
that the threshold flow rate needed to trigger the translo- 
cation of the chain is higher in wider channels. In or- 
der to quantify this behaviour, we define the threshold 
flow rate following the criterion P(J*) = P* = 0.5, 
and measure the corresponding values of J*, which, as 
shown in Fig. 02b), increases monotonically with D. In 
order to compare with the blob-model prediction, we 
carry out a fit of the simulation data to the function 



J*/J c = A/ln(X/D), obtaining A = 0.28 ± 0.01 and 
X = 6.88 ±0.08. 

From these fitting parameters and using Eq. ([5]) we 
obtain r m w 10 4 , which is consistent with the order of 
magnitude of the passage time of the chain through the 
nanochannel (« 5 x 10 4 timesteps). The predicted in- 
crease of the threshold flow rate with the channel width 
is thus confirmed by the numerical results, as shown in 
Fig. EJb). As expected, the data start to deviate from 
the prediction for D = 6 < R, where the chain begins to 
cross over to the unconfined regime, as found in the equi- 
librium simulations described above. It is thus pleasing 
that two very different approaches to coarse-graining, the 
blob and bead-spring models, agree for both the statics 
and dynamics presented here. 

This newly reported and unexpected dependence of the 
threshold flow rate for polymer translocation on the size 
of the nanochannel originates from the space-filling char- 
acter of the chains and their Brownian motion within the 
constriction. Given that the blobs composing the poly- 
mer match the width of the channel, the attempt fre- 
quency associated with the self-diffusion of the chain in 
the channel decreases with increasing pore width. This 
is both because larger blobs self-diffuse over a larger dis- 
tance, and because the viscous resistance to their motion 
scales with their size. It is counter intuitive that for fully 
confined polymers the probability of translocation is re- 
duced as the width of the channel is increased, leading 
to lower threshold flow rates in narrower channels. 

A feasible experimental setup to test our prediction, 
and to thereafter correctly interpret and understand ex- 
perimental data, follows from Beguin et al. [26j , who car- 
ried out measurements of the permeate flux and the re- 
jection coefficient (related to the probability of transloca- 
tion) for the injection of hydrosoluble polymers through 
nanopores etched in polycarbonate membranes. This 
corresponds to a blob size comparable to the size of the 
monomers, therefore making a direct comparison to ex- 
periments possible. A situation of great practical impor- 
tance is the confinement of single DNA chains in narrow 
channels, where the persistence length, Z p , is larger than 
the diameter of the confining channel, D. In the scaling 
description, the polymer length scale analogous to the 
blob size is the deflection length, l c ~ l^ 3 D 2 ^ 3, which, 
like the blob size, is an increasing function of D. Wc 
therefore expect that the same qualitative results for the 
attempt frequency hold in this situation. For larger pores 
(in the microfluidic regime) , where the persistence length 
of the chain is smaller than the size of the channel, the 
polymer behaves as a flexible chain, and is therefore also 
described by our theory. 

The threshold flow rate is a quantity of interest both 
from the biological and technological perspectives. Our 
findings show that a detailed description of the energy 
barrier crossing is still necessary, and we hope to inspire 
further research on potentially important effects arising, 
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for example, from polymer-pore interactions and pore ge- 
ometry. Our present results for the size dependence of 
the threshold flow rate may help in the rational design 
of filtering nanofluidic devices for which the channel size, 
rather than the operation flow rate, is engineered to trig- 
ger or hinder the translocation of the chains, and may 
shed light on possible mechanisms present in biological 
systems to control macromolecule transport. 
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